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£NJ ■ We analyze the phases of dense quark matter at nonzero quark 

5j . masses. We map out the phase diagrams for charged quark matter as 

02 [ well as charge neutral matter containing leptons. We find that in Color 

(— I ■ Flavor Locked (CFL) quark matter, the symmetric phase with equal 

*/""} ■ numbers of u, d and s quarks is unlikely to occur in nature. Various 
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less symmetric phases with additional meson condensates may play 



J> ' a role in neutron star cores, however, in agreement with the analysis 

^O . of Bedaque and Schafer. We discuss the implications of these novel 

^j color superconducting quark matter phases and their transitions for 

£""». ' neutron stars. 
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1 The meson condensed phases 

There has recently been renewed interest in the properties of dense quark 
matter |], [|. Three flavor QCD with massless quarks at high baryon den- 
sity is thought to have a particularly symmetric ground state, the Color 
Flavor Locked (CFL) phase [§]. In this phase the SU(3) co \ OT x SU(3) L x 
SU(3)r x U(1)b symmetry of QCD is broken down to the global diagonal 
SU(3) symmetry due to BCS-like pairing between quarks near the Fermi 
surface^. As a result, the gluons become massive via the Higgs mechanism. 
There is the possibility that color superconductivity is relevant to the struc- 
ture of neutron stars. However, in realistic situations, quark masses are 
nonzero and there may be nonzero chemical potentials for charge or lepton 
number, and the exact flavor symmetry of the QCD Lagrangian is reduced 
to U(l) em x U(l)y x U(1)b, in which case one might expect the CFL con- 
densate to preserve the U(l) em x U(l)y subgroup of SU(3) not explicitly 
broken?]. One can best understand how the symmetric CFL phase distorts 
under these applied stresses by studying the effective theory for the light 
excitations 0, || [7], §, ||, [K| |TTJ, the subject of this Letter. In particular, we 
explore the phase diagram for realistic quark masses in a density regime that 
may occur in neutron star interiors, exploiting recent observations by Schafer 
and Bedaque |12 |. We find that the symmetric CFL phase is unstable with 



respect to meson condensation, except above baryon densities approximately 
10 21 times that which might be expected in a neutron star core. Meson 
condensation results in a number of less symmetric color superconducting 
phases that may be relevant, with implications for the thermodynamic and 
transport properties of neutron star cores. 

The low lying spectrum of the symmetric CFL phase with massless quarks 
consists of a nonet of massless Goldstone bosons transforming under the 
unbroken SU(3) as an octet plus a singlet. A tenth light meson is also 
expected, a pseudo-Goldstone boson arising from spontaneously broken axial 

3 SU(3) co i or is the gauged symmetry of QCD; SU(3)l X SU(3)r is the global symmetry 
which acts on the three flavors of left- and right-handed quarks; U(\)b is the baryon 
number symmetry. 

4 Throughout this paper we will refer to both the SU(3) symmetric phase in the absence 
of quark masses and chemical potentials, and to the U{\) em x U(1)y symmetric phase in 
the presence of masses and chemical potentials, as the "symmetric CFL phase". The 
identifying characteristic of the symmetric CFL ground state is that it contains equal 
numbers of u, d and s H. 



U(1)a --an approximate symmetry broken by instanton effects, which are 
expected to be weak at high density. The massless mesons allow the system 
to respond continuously to applied stresses, such as the imposition of nonzero 
chemical potentials for SU(3) or U(1)b currents, since there is no gap in these 
channels for massless quarks. 

Turning on nonzero quark masses induces a gap in the spectrum. This 
would seem to make the symmetric CFL phase rigid against distortion due 
to applied chemical potentials. However, as recently pointed out by Schafer 
and Bedaque Jl2|, while quark masses produce gaps in the spectrum, they 



also induce new stresses on the system, acting in a manner analogous to an 
applied chemical potential and favoring meson condensation. As discussed 
by those authors, nonzero quark masses can be expected to lead to meson 
condensation and a complicated phase diagram. 

We can parameterize the effective theory describing excitations about 
the £77(3) symmetric CFL ground state in terms of the two fields B = 
e */3//s anc [ vj — qM^/U+v //a) ; representing the Goldstone bosons of broken 
baryon number (3, and of broken chiral symmetry, the pseudoscalar octet 7T, 
and the pseudo-Goldstone boson rj', arising from broken approximate U(1)a 
symmetry. Color degrees of freedom have been integrated out, and so under 
the original global SU(3)l X SU(3)r x U(1)b x U(1)a symmetry, these fields 
transform as the representations B = (1, 1) 10 and S = (3, 3)o,-4- 

The quark mass matrix M = diag(m u , m^, m s ) acts like a spurion, trans- 
forming as the representation (3, 3) 0j 2- Due to the U(1)a charges, oper- 
ators in the effective theory must involve even powers of M. Although 
there are two independent operators one can construct with two powers of 
M, the leading operator at weak coupling involves only the combination 
M = \M\M~ 1 = diag(mrfm s , m s m u , m u m<i), which transforms as a (3,3)o,4 
and can couple to X as TrMX. It is this operator that gives the leading 
contribution to the meson masses. However, one can also consider operators 
involving one power of M and one power of ML For example, the trace- 
less part of MLM, transforms as (8, l)o,o> the same as a left-handed flavor 
current. It has recently been pointed out by Bedaque and Schafer that the 
combination MLM/(2/x) (where \x is the quark number chemical potential) 
enters the effective theory exactly as would a chemical potential for flavor 
symmetries |12| f\ As they discuss, this term destabilizes the ground state 
for sufficiently large strange quark mass and leads to kaon condensation. 



5 A similar observation was made in Ref. fR 



The leading terms in the effective Lagrangian describing spatially con- 
stant nonet Goldstone boson fields ir and rj' in the presence of a chemical 
potential \xq for electric charge is given by 



£ = /, 



jTr D ED o tf + ]- f -A{D Q r}'f + ^Tr M (E + E+) + ^TrQEQEt 



4 " - 2/2^ "<> 2 " v / 2 

A)E = 3 E - ^ [0*qQ + X L ) E - E (// Q Q + X H )] . (1) 

The decay constants /„- and /^ have been computed previously 0. Q is 
the electric charge matrix diag(2/3, —1/3, —1/3) while X LR are the Bedaque- 
Schafer terms: X L = -^f- , X R = -^f , and M = \M\M~ 1 . 

The leading contributions to meson masses are the a and b operators. 
The coefficient a has been computed and is given by a = 3 -^-p [0 (f° r a 



discussion of possible log corrections see 0, [TO], [ll]]). The b term accounts 
for electromagnetic corrections to the charged meson masses. The coefficient 
b is proportional to a (the fine structure constant for the massless photon 
in the SU(3) symmetric CFL phase) and its value is scheme dependent^. 
From naive dimensional analysis, one can estimate the contribution to b from 
photon loops in the effective theory, leading to b ~ f- A 2 ; this formula reflects 
the quadratic divergence of a photon loop in the effective theory, which has a 
momentum cutoff equal to the gap A. This estimate is consistent with that 
of [p!4| , given that that reference did not attempt to estimate the factors of 



An . If the correct value of b differs from the crude estimate given above, our 
results will be changed quantitatively but not qualitatively. 
The meson masses in terms of the parameters a, b are 



ml- -- 


= a(m u + md)m s + b 


m\- -- 


= a(m u + m s )md + b 


m 2 K o - 


= a(m d + m s )m u . 



(2) 

The free energy density corresponding to the Lagrangian in eq. (|l|) in the 
mean field approximation, relative to the free energy density of the symmetric 
CFL phase (corresponding to E = 1) is: 

f 2 r 

n = J -f [Tr [/}, Ep, E f ] - 2 a Tr M(E + E f - 2) - b Tr [Q, E][Q, E f ] 



6 If one uses a mass independent subtraction scheme in Landau gauge, the one-loop 
electromagnetic corrections to the meson masses vanish, and the b term alone accounts for 
the entire electromagnetic splittings. We assume this scheme, and do not include explicit 
photon fields. 



(3) 

where S is a constant U(3) matrix characterizing the ground state. Working 
in the standard basis where the quark mass matrix is real, the matrix jl is 
given by 



jX = Hq + X , X = X L 



X, 



' 1\t 



(4) 



The effective chemical potentials vanish for the 7r°, rj and rj', while for the 
7r + , K + and K° mesons they are 
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The stationary points of the free energy Q with respect to variations of 
the meson fields are found as solutions to the matrix equation 



/ffit/iE -aME-b QE^QE 



— h.c. 







(6) 



To find solutions to this equation, we begin by considering E matrices cor- 
responding to a condensate of a single meson (0) = y/2 f n 9, where could 
be a charged pion, a charged kaon, or a neutral kaon field. Solving eq. (|JD 
yields stationary points with 9 ^ which for certain ranges of parameters 
have lower free energy than the symmetric CFL phase (6 = 0). We find the 
solutions 
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Evidently the nontrivial solutions with 9 ^ represent phases with lower free 
energy than the symmetric CFL phase. The charged pion and kaon phases 



(7) 



can have either electric charge, depending on the sign of appropriate effective 
chemical potential jl. We will refer to these 3-flavor, color superconducting, 
meson condensed phases as CFLn^, CFLK 1 * 1 and CFLK respectively, or 
simply as tc^, K and K°, distinguishing them from the more symmetric 
CFL phase. 

2 The phase diagram 

In Fig. 1 we plot the phase diagram for color superconducting quark matter, 
as a function of m 2 /(2/i), where /i is the quark number chemical potential, 
versus the chemical potential for positive electric charge, fiQ. In order to 
construct this diagram we used the quark mass ratios derived from chiral 
perturbation theory |T5|] , m u /m d = 1/2, m d /m s = 1/20, with the additional 
assumption m s = 150 MeV. For the parameter b in the chiral Lagrangian, 
we used the rough estimate b = |^A 2 . Finally, for the relation between the 
color superconducting gap A and /i, we used the perturbative expression 
from Ref.([|l6[]) down to /i = 10 GeV. At lower densities, we choose a smooth 
interpolant that mimics phenomenological models and predicts A ~ 100 MeV 
for /i ~ 500 MeV [|TJ. The superconducting gap obtained in this way is shown 
in Fig. |2|. We note that, even in the perturbative regime, the numerical value 
of the gap is still not well understood due to ambiguities associated with 
scale dependence of the gap equation |T7J . However, we find that although 
the precise location of the phase boundaries in Fig. [I] depend on the numerical 
value of the gap, the qualitative features of the phase diagram, which will be 
discussed below, are fairly independent of A(ju). 

Several features of the phase diagram Fig. 1 are striking. The first is that 
the symmetric CFL phase is disfavored until quark densities greater than 
~ 10 21 times what one might find in neutron star cores. Instead, at realistic 
values for /z, the K° phase is favored for small fiQ, while the K + and tt~ phases 
are favored for [Iq above or below critical values. The K° and K~ condensed 
phases are disfavored as they introduce strange quarks in a system trying 
to reduce their numberQ while for positive fiQ K + is favored over ir + both 
because the K + is lighter, and because a K + condensate reduces the strange 
quark number while a pion condensate does not. 



7 This is in contrast to kaon condensation in hadronic matter |l^, |l9|j , which suffers 
from a paucity of strange quarks, while the symmetric CFL phase has a surplus. 



As one approaches from either side the boundaries between the K° phase 
and the K + or it~ phases, one finds that the meson condensate angles 9 given 
in eq. ([7|) do not vanish. This indicates that the phase transitions are first 
order. One can verify this by solving the stationarity equation eq. (|6]) with 
a two condensate ansatz for the £ matrix. As an example we analyze the 
K + /K° transition. With the ansatz 
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one finds a saddle point solution to eq. 
K + phases, given by 
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This saddle point solution exists only for ranges of parameters for which the 
angles in eq. ([]) are real. The free energy density of this solution relative to 
that of the K° phase is 
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It is the quantity SQ which will control the tunneling rate between the K + 
and K° phases during a first order phase transition. 



3 Neutron stars 

If quark matter were to exist in neutron stars, the typical quark chemical 
potential is /i ~ 400 MeV. As discussed earlier, perturbative calculations are 
inapplicable here and model calculations based on phenomenological four- 
fermion interactions suggest that three flavor massless quark matter may be 
characterized by a gap as large as A ~ 100 MeV. As will become clear from 
the following discussion, this is larger than the typical energy associated with 



external stresses that quark matter is subject to in the neutron star interior. 
Therefore we can calculate the response of quark matter to these stresses 
within the frame work of the effective theory defined by eq. (1). 

Neutron stars are born in the aftermath of a core collapse supernova ex- 
plosion. The inner core of the newly born neutron star is characterized by 
a temperature T ~ 30 MeV and a lepton fraction (lepton number/baryon 
number) Y L ~ 0.3 implying /i e = /i„ e — fiQ ~ fi Ue ~ 200 MeV. The high tem- 
perature and finite lepton chemical potentials are a consequence of neutrino 
trapping during gravitational collapse p0| . Subsequently, the interior cools 



and looses lepton number via neutrino diffusion. This occurs on a time scale 
of few tens of seconds. At later times, T ^C 1 MeV, fi Ue = and the electron 
chemical potential is determined by the condition of electric charge neutral- 
ity. In noninteracting quark matter, fi e ~ m 2 s /2fi ~ 25 MeV. However, in the 
superconducting phase with A > ml/Afi our phase diagram Fig. 1 indicates 
that the ground state is the K° phase. The K° phase is electrically neutral 
without electrons. This result has far reaching consequences because, being 
devoid of electrons, the quark phase has no low lying charged excitations up 
to the mass of the lightest charged Goldstone boson. In the symmetric CFL 
phase, neutrality is guaranteed by the enforced equality of the u, d, and s 
quarks ||]. For the K° phase, charge neutrality is a fluke and is due to the 
fact that the K + is heavier than the K°. 

The core of a neutron star during the early stages of its evolution will be 
charge neutral, and will have nonzero lepton number. Assuming local elec- 
tric charge neutrality, the possible ground states among the phases we have 
discussed are K° condensation with all the leptons being neutrinos; K + con- 
densation with charge neutralizing electrons; or a mixed K + /K° condensate 
with electrons, corresponding to < (f) < n/2 in eq. (8). Our analysis at 
fixed jjlq found that the latter phase was not a minimum of the free energy, 
but a saddle point. Remarkably enough, when leptons are included and local 
charge neutrality is imposed, the saddle point solution of eq. (11) becomes 
the true minimum of the free energy over the range that it exists. Further, 
we find that it smoothly connects the K° phase with K + phase resulting in 
second order transition transitions. Note that eq. (0) remains valid and char- 
acterizes the stationary points of the free energy even when charge neutrality 
is imposed at fixed lepton chemical potential. 

A phase diagram depicting each of the three phases discussed above 
are shown in Fig. |3|. Electric charge neutrality is enforced in the phase 
containing the K + mesons by requiring that the electron number density 

8 



n e = —dQ, meson /djj,Q. This uniquely determines the electron chemical poten- 
tial and hence their contribution to the total free energy. In Fig. |3| we restrict 
the quark chemical potential to the range of relevance for neutron stars: 400- 
1000 MeV; and make the simplifying assumption that the gap A = 100 MeV 
over this range. 

From Fig. |3| we can infer that the ground state will contain electrons for 
fi Ue > 6 MeV . The superconducting phase admits electrons with ease since 
the K + mesons are not much heavier than the neutral K° mesons. We find 
that the electric charge chemical potential \xq is only a few MeV even when 
\i Ve ~ 200 MeV. Therefore the induced stress on the quark phase, character- 
ized by //k+ = /J-q + X, is always small compared to A. Consequently, unless 
excluded by high temperatures, quark matter in neutron star cores will be 
color superconducting even when leptons are trapped. As lepton number 
is lost via neutrino diffusion, we can expect a second order transition from 
the K + to the two condensate phase. Upon further deleptonization, when 
[i Ve ~ 5 — 10 MeV another second order transition from the K°/K + phase 
(containing electrons) to the K° phase (without electrons) is to be expected. 

Although it is premature to speculate precisely how this will influence 
various observable aspects of the early evolution of neutron stars we com- 
ment on a few important consequences of our findings. Firstly, since the 
only hadronic excitations in the CFL0 phases with energies below the gap 
are the pseudo Goldstone bosons with masses ~ 10 — 20 MeV the leptons are 
likely to dominate the transport processes and the matter specific heat. Dur- 
ing deleptonization we can expect dramatic changes to these fore mentioned 
properties. This is likely to affect various aspects of the early evolution of 
the neutron star including neutrino transport. 

If magnetic fields play a role during this early phase, we can expect field 
topology to change greatly during the early deleptonization phase since the 
initial K + and K°/K + phases are electrical superconductors they will expel 
the magnetic field in the core due to Meissner effect; while at later times when 
the system is characterized by a neutral K° phase and the U(l) em is restored, 
the Meissner affect disappears and the core now admits the magnetic field. 

Finally, since the neutrino free ground state contains no electrons, we note 
that the total lepton number emitted in neutrinos will be 20% — 30% larger 
than in conventional scenarios. This is particularly, exciting because current 



8 This depends on the value of b, for example, for b = e 2 A 2 (167r 2 times larger than 
the value we use) electrons appear for \x Ve > 20 MeV 



neutrino detectors such as SNO, which can distinguish between different fla- 
vors, in conjunction with larger detectors, such as Super Kamiokande which 
primarily detect the anti-electron neutrinos, is in principle sensitive to the 
lepton number emitted from a galactic supernova. 

It must be stressed that there are significant uncertainties in our analy- 
sis. Most importantly, the density regime one expects in neutron star cores is 
not high enough for the application of perturbative QCD. At these densities, 
the question of whether a quark phase exists at all, whether color supercon- 
ductivity occurs should such a quark phase exist, and what the size of the 
resulting gap and chiral Lagrangian parameters would be, are all presently 
unanswerable. In our analysis we have extrapolated from the perturbative 
regime and respected the symmetries of QCD, with the exception that we 
have assumed U(1)a violating instanton effects to be small. If the latter as- 
sumptions proves to be false at the relatively low densities found in neutron 
stars, a TrMS will be induced in the chiral Lagrangian || [K|, altering our 
analysis. We also note that we have ignored finite temperature effects in our 
analysis of the phase diagram. Thus we require that the matter temperature 
be small compared to the superconducting critical temperature T c ~ A . 
In addition, we have neglected finite temperature corrections to the effec- 
tive theory and the role of thermally excited mesons, which we expect to be 
small in the protoneutron star context. However, the contribution of thermal 
mesons to the free energy may not be negligible. Nevertheless, we feel con- 
fident that should color superconductivity exist in neutron stars, the phase 
diagram we have computed will be qualitatively correct and warrants further 
investigation. 
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Figure 1: Meson condensed phases in the neighborhood of the symmetric 
CFL state are shown in the (m|/2/x) — /Iq plane, where m s is the strange 
quark mass (set to 150 MeV), /i is the quark number chemical potential, 
and fiQ is the chemical potential for positive electric charge. At five times 
nuclear density \i ~ 400 MeV and (rn^/2/j,) ~ 25 MeV. Solid and dashed 
lines indicate first- and second-order transitions respectively. 
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Figure 2: A model for superconducting gap A as a function of the quark 
chemical potential /i. The curve agrees with perturbative calculations [T^ 
for fi > 10 GeV. 
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Figure 3: The phase diagram for charge neutral matter as a function of the 
quark number and neutrino number chemical potentials /i and \x Ve . The phase 
transitions shown are both second order. 
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